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ABSTRACT 



We have modelled the X-ray emission of T Tauri stars assuming that they have 
isothermal, magnetically-confined coronae. These coronae extend outwards until either 
the pressure of the hot coronal gas overcomes the magnetic field, or, if the corona 
interacts with a disk before this happens, by the action of the disk itself. This work is 
motivated by the results of the Chandra Orion Ultradeep Project (COUP) that show 
\ an increase in the X-ray emission measure with increasing stellar mass. We find that 

this variation (and its large scatter) result naturally from the variation in the sizes of 
the stellar coronae. The reduction in the magnitude of the X-ray emission due to the 
\& , presence of a disk stripping the outer parts of the stellar corona is most pronounced 

' for the lower mass stars. The higher mass stars with their greater surface gravities 

have coronae than typically do not extend out as far as the inner edge of the disk 
and so are less affected by it. For these stars, accretion takes place along open field 
I 1 lines that connect to the disk. By extrapolating surface magnetograms of young main 

sequence stars we have examined the effect on the X-ray emission of a realistic degree 
of field complexity. We find that the complex fields (which are more compact) give 
densities of some (2.5 — 0.6) x 10 10 cm~ 3 . This is consistent with density estimates of 
(1 — 8) x 10 10 cm~ 3 from modelling of individual flares. A simple dipole field in contrast 
gives densities typically an order of magnitude less. For the complex fields, we also 
find surface hotspots at a range of latitudes and longitudes with surface filling factors 
of only a few percent. We find that the dipolar fields give a relationship between X- 
ray emission measure and stellar mass that is somewhat steeper than observed, while 
the complex fields give a relation that is shallower than observed. This may suggest 
that T Tauri stars have coronal fields that are slightly more extended than their main 
sequence counterparts, but not as extended as a purely dipolar fields. 

Key words: stars: pre-main sequence - stars: activity - stars: imaging - stars: 
rotation - stars: spots 
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1 INTRODUCTION 

The recent Chandra Orion Ultradeep Project (COUP) has 
provided a remarkable dataset of the X-ray emission of pre- 
main sequence stars. It shows that for these T Tauri stars 
(as for main sequence field stars) the X-ray emission mea- 
sure rises slowly with stellar mass. For the T Tauri stars 
that show signs of active accretion, however, the correla- 
tion is less well defined UPreibisch et al. 20051 . These COUP 
stars have rotation rates (or Rossby numbers) that place 
them in the saturated or supersaturated part of the X-ray 
emission/rotation rate relation, but with emission measures 
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slightly below that of their main-sequence counterparts. 
This suppression of the X-ray emission also appears to be 
related to the presence of active accretion from the sur- 
rounding disk. The detection of significant rotational mod- 
ulation in some of the saturated stars suggests that the 
saturation of the X-ray emission is not the result of the 
complete coverage of the stellar surface in emitting regions 
JFlaccomio et al. 200511 . Rather, the emitting regions must 
be distributed unevenly around the star and they do not 
extend significantly more than a stellar radius above the 
surface. 

Observations of flares also yield information about the 
structure of T Tauri coronae. Although flares typical of the 
compact flares seen on main-sequence stars have been ob- 
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Figure 1. Magnetic structure of a dipole field that has been dis- 
torted by the outward pressure of the hot coronal gas trapped on 
the field lines. Beyond some radius (shown dotted) the pressure 
of the hot coronal gas is large enough to force open the mag- 
netic field lines. An example is shown of the last closed field line 
that intersects the equatorial plane at a radial distance r m and 
emerges from the stellar surface at a co-latitude m . If the coro- 
nal temperature and hence pressure is increased, or the stellar 
field strength is decreased, then r m moves inwards and the angle 
m increases. As a result, the fraction of the stellar surface that 
is covered in closed field capable of confining X-ray emitting gas 
decreases. 



served in the COUP sample, extremely large flares have also 
been detected. Estimates of plasma densities based on these 
flare observations are in the range of (1 — 8) x 10 10 cm -3 
( Fava ta et al. 20051 . The derived lengths of the flaring loops 
are around the location of the co-rotation radius ( 5R*) for 
the most reliable estimates, although some loop lengths ap- 
pear to much greater, suggesting perhaps that these flares 
have taken place in loops that link the star and the sur- 
rounding disk. 

The X-ray data from COUP therefore suggest that T 
Tauri stars have coronae that are structured on a variety 
of length scales, ranging from compact loops on scales 
of a stellar radius to much longer structures on the scale 
of the inner regions of an accretion disk. Observations 
of T Tauri magnetic fields also point to a very complex 
field structure, with localised regions of high field strength 
on the stellar surface. Unpolarized Zeeman-splitting 
measurements for photospheric lines with high Lande 
g-factors typically yield field strengths on T Tauri stars 
of order 2 to 3 kG (Johns- Krull, Valcnti & Koresko 1999 
|Johns-Krull, Valenti fc Saar*^0*0*4^ with surface fill- 
ing factors of order some tens of percent. Similarly, 
the observed splittings of the circularly-polarised Zee- 
man a components in the He I 5876 line indicate 



that the accretion streams impact the stellar surface 
in regions where the strength of the local magnetic 
field is of order 3 to 4.5 kG (IGuenther et al. 19991 

Johns- trull, Valenti fc Koresko 1999| 

Symington et al. 2005) . As |Valenti fc Johns-Krull (2004) | 
point out, unpolarised Zeeman-broadening measurements 
are sensitive to the distribution of magnetic field strengths 
over the stellar surface, but yield little spatial information. 
Circular polarization methods such as Zeeman-Doppler 
imaging give better spatial information about the magnetic 
polarity distribution, but have limited ability to determine 
field strengths at the level of the photosphere. 

The rate at which the field strength declines with ra- 
dial distance from the star depends on the complexity of 
the magnetic polarity distribution on the stellar surface. On 
the Sun, much of the flux emerging from small-scale bipolar 
regions in the photosphere connects locally to nearby re- 
gions of opposite polarity. The locally-connected field lines 
do not contribute to the coronal field at heights greater than 
the footpoint separation. As a result, the power-law depen- 
dence of the average coronal field on height steepens with in- 
creasing field complexity. A simple model assuming a global 
dipole with a surface field strength of order 2 to 3 kG would 
probably assume too high a value for the surface field, and 
allow it to decline too slowly with radial distance. 

Zeeman-Doppler images of post-T- Tauri stars give us a 
much better idea of the complexity of the magnetic field, 
and hence of the rate at which it falls off with height. The 
surface resolution of a Zeeman-Doppler image is determined 
by the ratio of the star's rotational line broadening to that 
of the intrinsic line profile convolved with the instrumen- 
tal resolution. On length scales smaller than a resolution 
element, the circular polarisation signals from small-scale 
field concentrations of opposite polarity cancel. The Zee- 
man broadenings seen in unpolarised light, however, do not. 
The "field strengths" in Zeeman-Doppler images are in effect 
magnetic flux densities, averaged over a surface resolution 
element (|Donati fc Collier Cameron 19971 . For this reason, 
the highest flux densities seen in Zeeman-Doppler images are 
always less than the maximum field strengths determined 
from unpolarised Zeeman splittings. Potential-field extrapo- 
lations from Zeeman-Doppler images should give a realistic 
measure of the field strength as a function of height beyond 
a few tenths of a stellar radius above the photosphere. 

The X-ray coronae of these pre-main sequence stars 
appear therefore to display a degree of spatial complex- 
ity and activity that is similar to that of their very active 
(saturated or supersaturated) main sequence counterparts, 
but with some significant differences. These differences (the 
suppression of X-ray emission and the greater scatter in 
the L x — Mi, relation) appear to be related to the pres- 
ence of active accretion. For main sequence stars the ob- 
served variation with rotation rate of both the magnitude 
of the X-ray emission its rotational modulation can be ex- 
plained by the centrifugal stripping of the outer parts of 
the corona lUnruh fc Jardine 19971 IJardine fc Unruh 1999 
I Jardine 2004t . This becomes important once the star is ro- 
tating sufficiently rapidly that its co-rotation radius is within 
the corona. |Ivanova fc Taam (2003) | have shown that this 
process can also reduce the angular momentum loss rate of 
rapid rotators (relative to more slowly rotating stars) with- 
out requiring a saturation of dynamo activity at high rota- 
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Figure 2. Contours of D (a measure of the goodness of fit to the 
COUP emission measures) as a function of the average value of 
the magnitude of the surface field strength and the constant of 
proportionality K which determines the pressure po a * the base of 
the corona since po = KBq. The cross marks the position of the 
minimum value of D and the lower panel shows in greater detail 
the area around this minimum. The plot shown is for models of 
the higher-temperature emission measures. A plot for the lower- 
temperature emission measures is qualitatively similar. 



tion rates. While most T Tauri stars are not rapidly rotating, 
their lower surface gravities and therefore greater pressure 
scale heights mean that their coronae may be more extended 
than their main-sequence counterparts. This suggests that 
a similar process may be important for these stars too. The 
presence of active accretion from the inner edge on an accre- 
tion disk may distort the corona of a T Tauri star, leading 
to a similar (but additional) stripping effect. 

The aim of this paper is to model the X-ray emis- 
sion from T Tauri stars, using the data (masses, radii, 
rotation periods, X-ray emission measures and tempera- 
tures) provided by the COUP results ijOetman et, al. 2 005a 
IGetman et, al. 2005bfl . We show these apparently puzzling 
features of X-ray emission from T Tauri stars can be ex- 
plained in the context of the varying sizes of stellar coronae 
that are limited either by the gas pressure of the stellar 
corona or by the presence of a disk. 



2 SIMPLE MODEL 
2.1 Passive disk 

We begin by addressing the question: what is the extent 
of the star's corona? We consider first the situation where 
each star has a disk that is passive in the sense that it does 
not distort the star's magnetic field. The coronal extent can 
then be estimated by calculating the height at which the gas 
pressure begins to exceed the magnetic pressure. In order to 
do this however we need a model for the star's magnetic 
field structure. As a simple example, we take a dipolar field 
geometry (see Fig. for which 
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B B = 



msinf 



(2) 



If we assume that the gas that is trapped on each closed 
field line is isothermal and in hydrostatic equilibrium then 

the gas pressure is p = p^e 1 ^ 1 9sds where g s = (g.B)/|B| 
is the component of gravity along the field. Here, for a star 
with rotation rate w, 



g(r,e)=(-GM i ,/r 
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The gas pressure at the footpoint of the field line po 
is a free parameter of this model. For this initial simple 
analysis we choose to scale this to the magnetic pressure of 
the average surface field strength Bo =< \B\ > such that 
po = KB where K is a constant that is the same on every 
field line. By scaling K up or down we can scale the overall 
level of the coronal gas pressure and hence the density and 
emission measure (see Appendix A for the calculation of the 
emission measure and the density). For each value of K we 
can also determine the variation of the gas pressure along 
each field line and hence the height beyond which the gas 
pressure exceeds the magnetic pressure. We show in Fig. 
the last closed field line, i.e. the most extended field line that 
can still contain coronal gas. Since the flux function 

A = ^ (4) 
r 

is constant along field lines, if this field line crosses the equa- 
tor at r = r m , then it has a footpoint on the stellar surface 
at a co-latitude O m given by 



sin 7i 



(•») 



If the field strength is decreased, or the temperature is in- 
creased, then the extent of the last closed field line will 
decrease and so m will increase. Our model now has six 
parameters: the mass, radius and rotation rate of the star, 
the coronal temperature and the values of the magnetic and 
gas pressures at the base of the corona. For a large num- 
ber of stars in the COUP sample the masses, radii and ro- 
tation rates have been determined. A two-temperature fit 
to the X-ray spectra then gives emission measures at two 
different temperatures. This leaves only the base gas and 
magnetic pressures undetermined. The value of po is sim- 
ply a scaling factor that raises or lowers the overall coro- 
nal pressure. We have chosen to make it proportional to 
the magnetic pressure since this gives a good agreement 
between predicted and observed values for emission mea- 
sures and mean densities both the the Sun and other main- 
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Figure 3. Emission measures for the high temperature (left) and low temperature (right) ranges from the COUP sample. Triangles show 
the observed values while crosses show the values calculated for a dipole field with values of the base magnetic and plasma pressures 
that give the best fit to the data. For each star we assume that any disk that may be present is passive in the sense that the maximum 
extent of the corona is limited only by the ability of the magnetic field to contain the hot coronal gas. The solid line is the best fit to the 
calculated values, while the dotted line is the best fit to the observed values, assuming the errors quoted for the log(emission measures). 
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Figure 4. Calculated emission-measure weighted densities for both the high temperature range (diamonds) and low temperature range 
(squares). In the left panel it has been assumed that any disks that may be present are passive so that the size of the corona is limited 
only by the ability of the magnetic field to contain the hot coronal gas. The right panel, in contrast, shows the densities obtained by 
assuming that the same stars each have instead an active disk so that the extent of the corona is limited by the action of the disk 
truncating the corona at the co-rotation radius. The lines show the best fits to the calculated densities. 



sequence stars (|Jardine, Collier Cameron fc Donati 2002| 
IJardine et al. 2002). We can then simply ask which combi- 
nation of (po, Bo) gives the best fit to the emission measures 
of the stars in the COUP sample. To do this, we select the 
116 stars from the COUP sample for which there are stel- 
lar radii, masses and rotation periods (fGetman et al. 2 005a 
IGetman et al. 2005bl . We exclude one star for which the 
apparent co-rotation radius is less than the stellar radius. 
For each star we then determine (for each of the two coro- 
nal temperatures quoted in the COUP sample) the emission 
measure (EM) and calculate, for the TV stars 

D = J? (logio(£ ; ^°b s .i)-log 10 (gM calc , i )) 2 (g) 
i=i a ' 

where the cr, are the errors in the log 10 (J5M) values. We 
show in Fig.|5]contours of D for a range of values of (po, Bq). 
Clearly, as we choose higher values of the base field strength 
(which will tend to give a larger corona) we need to have 



smaller values of the pressure (and hence the density) to 
give the same emission measure. Table^shows the values of 
(po,-Bo) that give the best fit to the data and Fig. |3]shows 
both the observed and calculated emission measure for these 
values. Even this simple model reproduces reasonably well 
the variation of emission measure with stellar mass and the 
degree of scatter at each value of M t . For every value of 
M* there are of course many possible combinations of the 
other parameters (i?*,P rot ,T) that determine the extent of 
the corona and the emission measure. If we fit a relationship 
of the form EM oc M* to the calculated emission measures 
we find that the slopes in both the high and low tempera- 
ture ranges are shallower than for a corresponding fit to the 
data, which are respectively 0.68 and 0.65 (assuming the er- 
rors quoted for the log(emission measures)). We remind the 
reader at this point that we are using only a subset of the 
stars from the COUP survey. We have restricted ourselves 
only to those that, in addition to emission measures and 
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Table 1. This shows, for a dipole field, the values of the aver- 
age field strength Bo =< \B\ > at the base of the corona and 
the corresponding gas pressure po = KB^ that give the best fit 
to the observed emission measures from the COUP sample. By 
fitting relationships of the form EM oc M% and n e oc Ml to the 
calculated values of the emission measures and densities we have 
determined the values of a and b. 
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Fig. shows the coronal extents (in units of the co- 
rotation radius for each star) that gave this best-fit value 
of (po,BZ). The more massive stars typically have coronae 
whose extents (relative to their co- rotation radii) are smaller 
than that of the lower mass stars. Indeed, for these field 
strengths and pressures many of the higher mass stars have 
coronae that do not extend as far as the co-rotation radius. It 
is the lower-mass stars with their lower surface gravity that 
are more likely to have coronae that extend out beyond the 
co-rotation radius and hence to suffer more interaction with 
any disks that may be present. This may have significant 
implications for the nature of the interaction between the 
stellar field and the disk and the exchange of torques be- 
tween them, since for these higher mass stars it may well be 
that the disk intercepts not the closed (dipolar) field of the 
star, but the open field of its wind. 



2.2 Active disk 
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Figure 5. Coronal extent (in units of the Keplerian co-rotation 
radius) as a function of stellar mass for stars in the COUP sample. 
Two stars (marked with a circle and square) have been selected 
as examples. Their coronal structure is shown more fully in Figs. 
E]and ED 



coronal temperatures, also have published masses, radii and 
rotation periods. For this subset of stars, the slope of the 
observed emission measure - mass relation is shallower than 
we would find if we were to use the full sample of emission 
measures. 

We can also (for each of the two temperatures quoted 
in the COUP database) calculate the emission measure- 
weighted density 



/ n 3 e dV 
JnldV 



(7) 



(see Fig. |1J. If we fit a relationship of the form n e oc 
to the calculated emission- measure weighted densities, we 
find that the slopes at the two temperature ranges are simi- 
lar, but the magnitude of the density for the lower tempera- 
ture range is systematically lower. As shown in the left-hand 
panel of Fig. these are rather lower than the values of 
(1 - 8) x 10 10 cm" 3 determined by |Favata et al. (2005)l from 
modelling flares observed in the COUP study, and indeed 
they are also lower than the values of 10 — 10 13 cm -3 cal- 
culated for young main sequence stars (for a comprehensive 
review see |Giidel . (2004^ . 



Up to this point, we have considered only the situation 
where each star has a passive disk. Each star's corona 
then extends out until the pressure of the coronal gas ex- 
ceeds that of the magnetic field confining it, causing the 
field to open up and release the gas into the stellar wind. 
For the lower-mass stars, however, whose coronae extend 
to distances greater than the co-rotation radius, the pres- 
ence of a disk extending in as far as the co-rotation ra- 
dius may disrupt this field. If this disruption is extreme, 
the presence of a disk may alter the structure of the stel- 
lar corona by shearing those field lines that pass through 
it, causing them to become open and therefore allowing 
coronal gas to escape in a wind |L ynden- Bell fc Boily 1994| 
|Lovelace, Romanova fc BisnovatyT^ogari^^95f ! We refer 
to this as an active disk. This process is likely to reduce 
the overall X-ray emission measure since it converts closed 
field lines that could contain hot X-ray emitting coronal gas 
into wind-bearing open field lines that will be dark in X- 
rays. This effectively places an upper limit on the size of 
the corona, since closed field lines cannot extend beyond 
the co-rotation radius. This would have the greatest impact 
on the lower-mass stars whose coronae are naturally more 
extended because of their lower surface gravity. The X-ray 
emission from stars whose coronae did not extend as far as 
the co-rotation radius would be affected much less severely. 

We therefore take the extent of the corona to have a 
maximum value of the location of the inner edge of the disk 
which we take to be at the co-rotation radius We then cal- 
culate the emission measure of the gas contained in each 
stellar corona for a range of values of both the mean mag- 
netic field at the stellar surface and the base pressure of the 
corona. We show in Table Q the values of the base gas and 
magnetic pressure that give the minimum value of D (as 
defined by ©) and show in Fig. |S|both the observed and 
calculated emission measures at these best-fit values. The 
presence of an active disk clearly gives a rather better fit 
to the emission measures at both temperature ranges. This 
is mainly because of the lower mass stars which have coro- 
nae that extend beyond the co-rotation radius and whose 
emission measure is suppressed by the action of the disk 
stripping the corona. As a result, these stars have a lower 
emission measure than would be produced with a passive 
disk and so the slope of the fitted line is larger and closer 
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Figure 6. Emission measures for the high temperature (left) and low temperature (right) ranges from the COUP sample. Triangles show 
the observed values while crosses show the values calculated for a dipole field with values of the base magnetic and plasma pressures 
that give the best fit to the data. We assume that each star has an active disk in the sense that the maximum extent of the corona is 
limited by the inner edge of the disk which is assumed to be at the co-rotation radius. The solid line is the best fit to the calculated 
values, while the dotted line is the best fit to the observed values, assuming the errors quoted for the log(emission measures). 



to that found for the observed emission measures. While the 
calculated densities (as shown in Fig. |1J are still lower than 
observed, they are a little higher than if the disk is assumed 
to be passive and the slope of their variation with mass is 
greater. This is because the best-fit value of the base pres- 
sure (and hence the base density) is higher for the an active 
disk where the typically smaller sizes of the stellar coro- 
nae require higher densities to match the observed emission 
measures. The lower-temperature coronae with the smaller 
scale heights are less affected by the presence of the disk and 
their density typically increases by less when an active disk 
is present. 



2.3 Surface hot spots 

Observations of magnetic fields on T Tauri stars and 
modelling of accretion signatures suggest that the frac- 
tion of the surface area of the star that is covered in 
accreting hot spots is very small, of the order of a 
few percent Ipuzerolle et al. 20031 ICaTvet et al. 20041 
IMuzerolle et al. 20051 IValenti fc Johns-Kmll "20041 

|Symington et al. 2005^ . We can place an upper limit 
on this by determining which field lines (in a corona whose 
size is limited only by pressure balance) would intersect 
the disk. Some subset of these field lines will be capable of 
supporting an accretion flow. We can determine this surface 
area fraction F by calculating the area of the two annuli 
on the stellar surface that lie between O m (the footpoint of 
the last closed field line) and Qkcr (the footpoint of the 
field line that passes through the co-rotation radius). As a 
fraction of the total stellar surface area this gives 



F = cos 6 m — cos 6kcr- 



(8) 



This fraction is shown in Fig.|7|as a function of the coronal 
extent. If the corona does not extend out as far as the inner 
edge of the disk, then only those open field lines that pass 
close to the equator will intersect the disk. Since these are 
likely to come from the polar regions, their area coverage 
at the stellar surface is very small. For those coronae that 
do extend beyond the inner edge of the disk, the area of 
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Figure 7. Crosses show the upper limit to the fraction of the 
stellar surface area that is covered in field lines that pass through 
the disk. A dipole field has been assumed. 



the surface covered by those field lines that cross the disk 
depends on both the value of the co-rotation radius and also 
the maximum extent of the corona as determined by pressure 
balance. This can be seen from J3J since, if the stellar radius 
is sufficiently small compared to either the co-rotation radius 
or the maximum extent (r m ) of the corona, 



F 



r KCR 



(9) 



This suggests that it is the lower-mass stars with their more 
extended coronae that should show the largest hotspots. 

It seems unlikely however that T Tauri stars will have 
simple aligned dipolar fields. The observations of signif- 
icant rotational modulation of the X-ray emission from 
many of the stars in the COUP sample are a clear in- 
dication that their fields are more complex than a sim- 
ple dipole. The predominance of a period for the X-ray 
emission that is equal to (or half) that of the photomet- 
ric period suggests the presence of either one dominant re- 
gion of X-ray emission on the stellar surface or two at op- 
posite longitudes. This is very similar to the predictions 
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Figure 8. Calculated coronal structure based on an example surface magnetic field distribution for one of the lower mass stars (shown 
as a circle in Fig.|3) where it is the presence of the disk that is primarily limiting the extent of the corona. In this example the stellar 
mass is 0.15 Mq, the radius is 4.02 Rq, the rotation period is 17.91 days and the natural extent of the corona (in the absence of a disk) 
would be at most 1.58 times the co-rotation radius. The left panel shows the field lines that could support an accretion flow, drawn from 
the inner edge of the accretion disk which is situated at the co-rotation radius. The middle panel shows examples of the field lines that 
are closed and contain X-ray emitting gas and the right panel shows the corresponding structure of the X-ray corona. 




Figure 9. Calculated coronal structure based on an example surface magnetic field distribution for one of the higher mass stars (shown 
as a square in Fig. [5]) where it is the pressure of the hot coronal gas that is primarily limiting the extent of the corona. In this example 
the stellar mass is 2.57 Mq the radius is 1.74 Rq, the rotation period is 1.38 days and the natural extent of the corona (in the absence 
of a disk) would be at most only 0.26 times the co-rotation radius. The left panel shows the field lines that could support an accretion 
flow, drawn from the inner edge of the accretion disk which is situated at the co-rotation radius. The middle panel shows examples of 
the field lines that are closed and contain X-ray emitting gas and the right panel shows the corresponding structure of the X-ray corona. 



of the X-ray rotational modulation for AB Dor based on 
extrapolation of surface fields determined from Zeeman- 
Doppler imaging ( Jardine, Collie r Cameron & Donati 2002 
I Jardine et al. 2002t . These extrapolations show that the 
large-scale field of AB Dor resembles a tilted dipole, giving 
(in addition to many small-scale bright regions) two domi- 
nant bright regions at opposite longitudes. The latitudes of 
these bright regions can be related (crudely) to the tilt of 
this large-scale dipole. Depending on the stellar inclination 
one or both may be visible. We therefore explore the impact 
on the X-ray emission from these T Tauri stars that a field 
with a realistic degree of complexity might have. 



3 FIELDS WITH A REALISTIC DEGREE OF 
COMPLEXITY 

Although the observations to date suggest that T Tauri 
stars have fields that are more complex than a simple 
dipole, we do not as yet have any direct observations of 
the detailed geometry of the field. We therefore use fields 
extrapolated from surface magnetograms of two young 
solar-like stars, AB Dor (period 0.514 days) and LQ Hya 
(period 1.6 days) that have been obtained using Zeeman- 
Doppler imaging (Donati et al. 1997 IDonati et al. 20031 
IDonati fc Collier Cameron 19971 IDonati et al. 19991 . This 
allows us to determine the influence of field complexity 
and to assess the effect that this might have on our results 
based on simple dipolar fields. Since the rotation axes of 
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both of these stars are inclined to the observer, only one 
hemisphere can be imaged reliably. In order to provide a 
realistic surface map for the hidden hemisphere, we use 
a map taken from another year. Thus for the first map 
we use magnetograms of AB Dor generated from data 
acquired in December 1995 and December 1996, and for 
the second map we use magnetograms of LQ Hya generated 
from data acquired in December 1993 and December 
2001. The coronal fields are extrapolated from these 
using the "Potential Field Source Surface" method orig- 
inally developed by |Altschuler fc Newkirk, Jr. (1969)] 
for extrapolating the Sun's coronal field from so- 
lar magnetograms. We use a code originally devel- 
oped by |van B allcgooijcn, Cartledge fc Priest (199*8)1 
Since the method has been described in 
|Jardine, Collier Cameron fc Donati (2002)] we provide 
only an outline here. Briefly, we write the magnetic field 
B in terms of a flux function <3/ such that B — — V\& 
and the condition that the field is potential (V x B = 0) 
is satisfied automatically. The condition that the field is 
divergence-free then reduces to Laplace's equation V 2 ^ = 
with solution in spherical co-ordinates (r, 6, <f>) 

N I 

* = J2 X! ^ aimT ' + b ^ Ml+1) ] P ^(0)e m ^ (10) 

1 = 1 m= — l 

where the associated Legendre functions are denoted by P; m . 
The coefficients a; m and bi m are determined by imposing the 
radial field at the surface from the Zeeman-Doppler maps 
and by assuming that at some height R„ above the surface 
(known as the source surface) the field becomes radial and 
hence Be(R s ) = 0. This second condition models the effect 
of the plasma pressure in the corona pulling open field lines 
to form a stellar wind. 

For rapidly-rotating solar-like stars, we can use the ob- 
served locations of the large slingshot prominences that typ- 
ically form in their coronae to estimate the radii at which 
the corona must still be closed. For T Tauri stars, however, 
we do not have this option. As a fairly conservative estimate 
we choose to take the radial distance at which a dipolar field 
(of the same average field strength) would remain closed. 
We choose the dipolar field as it gives the greatest coronal 
extent (since of all the possible multipolar components, the 
dipolar falls off most slowly with radial distance). If this 
maximum extent exceeds the co-rotation radius, we set it to 
be the co-rotation radius instead. 

We determine the pressure at every point by calculat- 
ing the path of the field line through that point and solving 
hydrostatic equilibrium along that path. Once again, we set 
the gas pressure at the base of the field line to be propor- 
tional to the magnetic pressure there, noting however that 
now the base pressure varies across the stellar surface. By 
varying the constant of proportionality K we can raise or 
lower the overall level of the gas pressure. If along any field 
line we find that the gas pressure exceeds the magnetic pres- 
sure, we assume that the field line should have been opened 
up by the gas pressure and so we set p — 0. 

We assume that the disk occupies a range of latitudes 
extending some 10° above and below the equator. We do not 
set the location of the inner edge of the disk, but allow it to 
adjust locally to the structure of the magnetic field at each 
latitude and longitude. We therefore assume that the disk 



Table 2. This shows, for two example field structures with a 
realistic degree of complexity, the average field strength Bq =< 
|S| > at the base of the corona and the corresponding gas pressure 
po = KBq. By fitting relationships of the form EM oc M° and 
n e tx M\ to the calculated values of the emission measures and 
densities we have determined the values of a and b. 



AB Dor-like LQ Hya-like 

T High Low High Low 

< \B\ > [G] 143 143 143 143 

po [dyne cm' 2 ] 3.6 0.2 13.2 0.9 

a 0.55 0.53 0.54 0.53 

b -0.05 -0.11 -0.06 -0.12 



extends in towards the surface of the star until it reaches the 
first field line along which a) the gas pressure is less than the 
magnetic pressure everywhere along the field line and b) the 
net gravitational force at the point where the field line passes 
through the disk acts inward towards the star. We assume 
that this first field line carries all the accreting material from 
that particular latitude and longitude and that it shields any 
field lines interior to it. Figs. [5] and [5] show the structure of 
these accreting field lines, examples of the closed field lines 
that could contain hot coronal gas and X-ray images from 
the same perspective. 

Of our two free parameters, one (the flux density or 
field strength) is determined by the surface magnetograms. 
As before, the other parameter is the constant of propor- 
tionality K that determines the gas pressure at the base of 
each field line. For the dipole field, we chose for simplicity 
to have a uniform po based on the average field strength 
< |B| >. For these complex fields, however, the base mag- 
netic pressure can vary significantly across the surface and 
as a result, so can the base gas pressure. In TableHtherefore 
it should be noted that the value of po = K < \B\ > 2 is just 
the value based on the average value of < \B\ >. 

We can therefore, for each of the two example field 
structures, determine which value of K gives the best fit 
to the emission measures from the COUP survey (see Fig. 
11011 . We show in Fig. llll the calculated emission measures su- 
perimposed on the observed emission measures. As can be 
seen from Table [2] the slopes for these calculated values are 
shallower than for the observed emission measures, but the 
values for the two field structures are remarkably similar. 

There is more of a variation to be found in the densi- 
ties. For the dipole field, densities of around 0.4 x 10 10 cm~ 3 
are typical for the high temperature range. For these more 
complex fields, however, the densities are almost an order of 
magnitude higher (2.5 x 10 10 cm -3 ). The slopes are also less 
than for the active disk with a dipolar field. This is probably 
because the more compact fields can support a higher den- 
sity, but they do not extend out as far as the dipolar fields 
and so they are less affected by the presence of a disk. 

With these more complex fields we would expect to find 
hotspots that are in different surface locations to those for 
a dipole field and that are perhaps different in size. Fig. 1131 
shows the fraction of the stellar surface area that is covered 
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Figure 10. Variation of D (a measure of the goodness of fit to the data) with the base pressure for (left) the field structure shown in 
Fig. |H| and (right) the field structure shown in Fig. |5| The solid lines denote the higher-temperature coronae, and the dashed line the 
lower temperature coronae. 





Figure 11. Observed emission measures for low (left) and high (right) temperature ranges from the COUP sample. The top panels refer 
to the field structure shown in Fig. [S] while the bottom panels refer to the field structure shown in Fig. [5] The solid line is the best fit to 
the calculated emission measures, while the dotted line is the best fit to the observed emission measures. 



in field lines that pass through the inner edge of the disk. 
For the higher mass stars whose coronae do not reach out 
to the disk, the field lines that reach the disk are open and 
radial. Variations in the position of the co-rotation radius 
then have little effect in determining which field lines pass 
through the disk and so the fraction of accreting field lines 
varies little with stellar mass. For those lower-mass stars 
whose coronae extend out beyond the co-rotation radius, 



the stellar magnetic field is still closed at the inner edge of 
the disk. Consequently, changes in the co-rotation radius 
can make large changes to whether a field line is accreting 
or not and so the surface area of hot spots shows a greater 
scatter. In either instance, however, the filling factor is 
low, in agreement with results from modelling of accretion 
signatures (IMuzerolle et al. 20031 ICalvet et al. 20041 
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Figure 12. Emission-measure weighted densities for both the high temperature (diamonds) and low temperature (squares) ranges from 
the COUP sample. The left panel refers to the magnetic fields shown in Fig, 151 and the right panel to the magnetic fields shown in Fig, 151 
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Figure 13. Surface filling factor of hotspots for (left) the magnetic fields shown in Fig. Inland (right) the magnetic fields shown in Fig. 151 



IMuzerolle et al. 20051 IValenti fc Johns-Krull 2004 

|Symington et al. 2005} . 



4 SUMMARY AND DISCUSSION 

We have shown in this paper that it is possible for a simple 
model to reproduce the observed variation with mass of the 
X-ray emission measure of T Tauri stars. This variation is 
due to the restriction in the size of the stellar corona that 
results either from the pressure of the hot coronal gas in 
opening up the magnetic field lines, or from the shearing 
effect of a circumstellar disk. The ratio of the pressure scale 
height to the co-rotation radius (a reasonable estimate of 
the location of the inner edge of an accretion disk) scales as 
A/i?KCR oc Mr 2 for a polytropic star. Hence, lower mass 
stars are more likely to have coronae that reach the inner 
edge of any accretion disk present and so their coronal ex- 
tent will be set by the co-rotation radius. Higher mass stars 
have coronae that may not extend as far as the co-rotation 
radius and so their coronal extent is set simply by the ability 
of the magnetic field to contain the gas. For these stars, ac- 
cretion will occur principally along open, nearly radial mag- 
netic field lines. 

Even in the absence of a disk, the variation of the ex- 



tent of the corona will produce a drop in the X-ray emis- 
sion measure with decreasing stellar mass. This process 
of coronal stripping has already been shown to reproduce 
the observed variation with rotation rate of both the mag- 
nitude and rotational modulation of the X-ray emission 
measures of main-sequence stars llUnruh fc Jardine '19971 
IJardine fc Unruh 19991 I Jardine 20041 . The stars in the 
COUP sample appear to lie in the saturated or super- 
saturated part of the rotation-activity relation (i.e. where 
X-ray emission measure is independent of rotation rate 
or even falling with increasing rotation rate). Their emis- 
sion measures are slightly lower however than their main- 
sequence counterparts and they show a greater scatter. 
JPreibisch et al. 2005t have examined this question of the 
overall suppression of the X-ray emission. They divide the 
sample of stars into those that show some indication of the 
presence of a disk (i.e. an infra-red excess) but no apparent 
signatures of active accretion (as measured by the width of 
the H a line) and those that show indications of both. The 
main contributors to the drop in X-ray emission appear to 
be those stars that show signs of active accretion and for 
which there is therefore clear evidence that their disks are 
interacting with the stellar magnetic field. Stars lacking this 
clear signature of disk interaction appear to show a level of 
X-ray emission that is similar to their saturated or supersat- 
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urated main sequence counterparts. This is consistent with 
our model which predicts that all stars will show a drop 
in emission with decreasing mass due to the pressure strip- 
ping of their coronae, but that in those stars where a disk 
is actively interacting with and distorting the stellar mag- 
netic field, this drop in emission will be further enhanced. 
|Preibisch et al. (2005) | also note that the large scatter in 
X-ray emission is most pronounced in the actively accret- 
ing stars. This can be understood in terms of our model 
since in those stars where the disk is actively stripping the 
corona, the coronal extent depends on a further parameter, 
the rotation rate, since this determines the location of the 
co-rotation radius where the inner edge of the disk is most 
likely to be situated. 

The observations of the strength of the magnetic 
field at the base of the accretion funnels and across 
the rest of the stellar surface (Valenti & Johns-Krull 2004 
|Symington et al. 2005^ combined with the observation of 
significant X-ray rotational modulation in stars in the 
COUP sample ( Flaccomi o et al. 200511 have prompted us to 
investigate the importance of the nature of the stellar mag- 
netic field. We have modelled the extent of the stellar corona, 
the X-ray emission measure and the density both for a sim- 
ple dipole field and also for fields with a realistic degree of 
complexity. We do not yet know in detail the structure of 
T Tauri coronal fields, but we have used instead coronal 
magnetic fields extrapolated from surface magnetograms of 
young main-sequence stars. While we do not yet know if 
these capture the coronal structure of T Tauri stars, they 
allow us to explore the effect of a more compact field. This 
is important as dipole fields (the most commonly used ex- 
ample until now) fall off with height more slowly than any 
other type of field except a purely radial one. They therefore 
model the most extended type of field possible. Fields that 
are more complex are more compact: their field strengths 
decay with height more rapidly than a dipolar field and so 
they will typically be less extended than a dipole and con- 
sequently are likely to produce a higher density. 

We find that the more complex fields do indeed produce 
more compact, higher density coronae, with densities in the 
range (2.5 — 0.6) x 10 10 cm -3 for coronae whose temperatures 
are set to the "high" and "low" temperatures in the COUP 
database. These densities are nearly an order of magnitude 
higher than for a dipole field. The densities produced by 
different types of complex fields are relatively similar, with 
the main differences being the the locations of the surface 
hotspots (where the accretion flows impact on the stellar 
surface) and their filling factors. In contrast to dipole fields 
where any accretion flow must reach the surface at the 
poles (thereby producing no rotational modulation), more 
complex fields can produce hotspots at lower latitudes. 
|Mahdavi fe Kenyon (1998)| have modelled the locations of 
hotspots in dipole field that are tilted with respect to the 
rotation axis, but these can only provide one hotspot in each 
hemisphere, while more complex fields can produce many 
smaller hotspots at a range of latitudes and longitudes. We 
find filling factors are small (typically a few percent) in 
agreement with predictions made from modelling accretion 
signatures l|Muzerolle et al. 20031 ICalvet et al. 20041 
IMuzerolle et al. 20051 I Valenti fc Johns-Krull "20041 

|Symington et al. 2005) . 

We have also calculated the form of the relationship 



EM oc M+ based on the published stellar masses, radii, ro- 
tation rates and coronal temperatures for the COUP stars 
ICetman et al. 2005al IGetman et al. 2005bt . For the higher 
temperature range in the COUP sample, the observed emis- 
sion measures (allowing for the published errors) give a value 
of a — 0.68. If we assume that all stars have either no disk 
present, or only a passive disk that does not distort the stel- 
lar corona, we find a value of a = 0.34 for a dipole field, 
compared with a = 0.70 if we assume that all stars have an 
active disk that distorts any stellar corona that extends as 
far as the inner edge of the disk (which is assumed to be at 
the co-rotation radius). If we assume a more complex field 
geometry this value drops to a — 0.55 if the stellar disks are 
assumed to be active. A similar result is found for the low 
temperature component. Thus the presence of disks capable 
of distorting those stellar coronae that extend beyond the co- 
rotation radius seems to be necessary to explain the slope 
of the observed variation with mass of the X-ray emission 
measure. The degree of complexity of the field has a smaller 
effect on the overall X-ray emission, but it does have a sig- 
nificant effect on the densities. Assuming that all stars have 
dipole fields gives densities much lower than those derived 
from observations of flares j I av al a et al. 20051 . while the 
assumption of a more complex field geometry gives densi- 
ties comparable to the published values. These results taken 
together suggest that the true T Tauri stellar fields may 
be more extended than those of the main sequence counter- 
parts, but not by as much as a purely dipole field. They may 
simply have a dipolar component that is somewhat stronger 
than in the main sequence stars. 

In summary, we have shown that a simple model (which 
reproduces the X-ray activity- rotation relation for main se- 
quence stars) can explain the variation of X-ray emission 
measure with stellar mass observed in the COUP survey, 
if we allow for the presence of an accretion disk which will 
limit the size of any stellar corona that reaches out to its 
inner edge. 
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6 APPENDIX 

In deriving expressions for the variation of the emission mea- 
sure and density with radial distance from the star, we begin 
with . If we scale all distances to a stellar radius, we can 
write 



p(r,0)=p*(l,0)fflcp 



- l)+$e sin 2 8(r 2 -1) 



(11) 



where p*(l,#) is the pressure distribution across the stellar 
surface and 



GM/n 
9 k B T/m 



m 2 r 2 /2 
k B T/m 



(12) 
(13) 



are the surface ratios of gravitational and centrifugal ener- 
gies to the thermal energy. The emission measure Em — 
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J n^dv is then 



EM = 



47rr 2 
(fc B T)2 



(14) 



/ p 2 (l,0)exp [-2$ c cos 2 <9(r 2 - I)sin0d0] 

Je=e m 

I 



9=e 

r = r m 



r 2 exp 



2$ s (--l + 2$ c (r 2 -l) 



If we now use the substitution fx = cos 9 such that ii m = 
cosO m = (1 — r/rm) 1 ^ 2 then the ^-integral can be written 

as 

f 1 p 2 (l,0)exp[-2<£> c/ x 2 (r 2 -l)]d M (15) 
Jo 

which, on using the substitution t = [2"l> c (r 2 — l)] 1//2 /i gives 



2 $ c ( r 2_ 1)1/2 J o 



pl(l,9)e- t2 dt 



(16) 



where t m = [2$ c (r 2 - 1)] 1/2 [1 — r/r m ] 1/2 . If we now assume 
that the pressure is uniform over the stellar surface we ob- 
tain the following expression for the emission measure as a 
function of r: 



EM(r) = 




'r=l 

exp 



2<£> c (r 2 -l)(l--^) 

' 5 



(17) 

1 1/2 



2$ 9 (i-l)2<I> c (r 2 -l) 



dr. 
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We can, in a similar way, obtain an expression for the emis- 
sion measure weighted density n e (r) = J n z e dV/ J n 2 e dV . 



n e (r) = 



EM(r) 

„ r _ r 




7T 3 2r'ipi 



3$ c (k B T)3 



l r=1 (^ 2 - 1) 1/2 
exp 



erf 



3$ c (r 2 - 1)(1 



(18) 

1/2 



3<& 9 (i-l)3<E. c (r 2 -l) 



dr 



where we remind the reader that all distances r have been 
scaled to a stellar radius. 



REFERENCES 

Altschulcr M. D., Ncwkirk, Jr. G., 1969, Solar Phys., 9, 131 
Calvet N., Muzerolle J., Bricefio C, Hernandez J., Hartmann L., 

Sauccdo J. L., Gordon K. D., 2004, aj, 128, 1294 
Donati J.-F., Collier Cameron A., 1997, MNRAS, 291, 1 
Donati J.-F., Semel M., Carter B. D., Rees D. E., Collier 

Cameron A., 1997, MNRAS, 291, 658 
Donati J.-F., Collier Cameron A., Hussain G., Semel M., 1999, 

MNRAS, 302, 437 
Donati J.-F. et al., 2003, MNRAS, 345, 1145 
Favata F., Flaccomio E., Reale F., Micela G., Sciortino S., 

Shang H., Stassun K. G., Feigelson E. D., 2005, apjs, 160, 

469 

Flaccomio E., Micela G., Sciortino S., Feigelson E. D., 

Herbst W., Favata F., Harndcn F. R., Vrtilck S. D., 2005, 
apjs, 160, 450 



